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errors  in  inspection  and  grading  of  samples  from  finite,  possibly 
stratified  lots  are  obtained.  Screening,  and  hicrarchal  screening 
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on  the  advantages  of  these  techniques  assessed. 


Accession  For 

NTIS  GRA&I 

DTIC  TAB 
Uii.0.r>nouncrd  □ 

Justificnti  or. _ 


f 


Thia  document  lv  -j  t'-^a  iryo  roved  ^ 
lor  public  roltKie ,  v  J  bulu;  iia 

dtottlbutton  la  unlimited. 


i 


1 .  INTRODUCTION 
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For  convenience,  we  first  summarize  some  results  obtained  by  Johnson 
et^  al_.  (1980)  and  Johnson  5  Kotz  (1981). 

A  random  sample  of  size  n  is  chosen  (without  replacement)  from  a  lot 
of  size  N,  which  contains  0  defective  (or  nonconforming)  members.  Suppose 
that  on  inspection  of  items  in  the  sample,  the  probability  that  a  defective 
item  will  be  classified  as  defective  is  p,  while  the  probability  that  a  non¬ 
defective  item  will  (erroneously)  be  classified  as  "defective"  is  p'. 

If  Y  denotes  the  actual  unknown  number  of  defective  items  in  the  sample, 
the  distribution  of  Y  is  hypergeometric  with  parameters  (n,  D,  N) .  If  X 
denote  the  number  (among  these  Y)  which  are  (correctly)  classified  as 
defective,  and  X'  the  number  (among  the  fn-Y)  nondefective  items  in  the 
sample)  which  are  (incorrectlyl  classified  as  "defective"  then,  conditional ly 
on  Y,  the  variables  X  and  X'  are  independent  binomial  with  parameters  (V,p'v 
and  (n-Y,  p’)  respectively.  Averaging  over  the  distribution  of  Y,  the 
distribution  of 

:  =  X  + X’ , 

the  total  number  of  items  described  as  "defective"  after  inspection  of  the 
sample,  can  be  formally  expressed  as 

Binomial  (Y,p)  *  Binomial  (n-Y,  p’)  y  Hypergeometric  (n.ll.N)  (Jj 

The  symbol  *  stands  for  convolution,  and  the  symbol  y  indicates  the  "correspond¬ 
ing"  operation  with  Y  distributed  as  "(srf.  e.g.  Johnson  f,  kot:  (I960, 

Chapter  8)).  Distribution  (1)  is  a  mixture  of  convolutions  of  two  binomial 
distributions.  The  r-tli  descending  factorial  moment  of  Z  is 
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where  Z^  =  Z(Z-l)  ...  (Z-r+1).  In  particular 
E(Z)  =  np/N 

—r,  n(n-l)  D  ,,  nl ,  .. 

var(Z)  =  np(i-p)  -  — y2-  *  ^  O-d(p-p’) 


(3.1) 

(3.2) 


where  p  =  {Dp  +  (N-D)p’)N  is  the  probability  that  a  single  individual, 
chosen  at  random  from  the  lot,  will  be  described  as  "defective"  after 


inspection  (whether  it  is  really  defective  or  not).  Tables  of  the 


distribution  of  Z  for  n=10  with  p  =  0.75(0.05)0.95;  p'  =  0(0.025)0.1; 

N  =  100  and  D  =  5,  10,  20;  N  =  20o  and  !)  =  10,  20,  10  are  presented  in 

Johnson  and  Kotc  (1981).  If  p'  =  0,  wc  have  the  situation  described  in 

Johnson  et  al_.  (1980)  and  (1)  becomes  a  hypergeometric-binomial  distribution. 

For  fixed  values  of  D/N  =  X  say,  the  distribution  (1)  is  quite  sensitive  to 

changes  in  p  and  p',  but  not  to  changes  in  1)  and  N  unless  n/N  (the  sampling 

fraction)  is  large.  It  is  ea.-y  to  see  that  as  D,  N  ->  «  (with  D/N  =  X j  the 

distribution  of  Z  tends  to  binomial  with  parameters  (n,Xp+l (l-X)p’) . 

Johnson  and  Kotc  (1981)  also,  alia,  consider  lots  divided  into  k 

strata  IT  ,  . ...  U.  of  sices  N  ,  N.  (with  )  N.  =  N)  and  suppose  that  the 

i  K  1  k  .  3 

i=i 

probability  of  an  individual  from  the  j-th  stratum  being  classified  as 


defective  is  p^ .  (The  situation  considered  at  the  beginning  of  this  section 
corresponds  to  k=2,  N=D,  N -,=N - H ,  p,  =  ]>,  p,  =  p'  .)  The  distribution  of  the 
total  number  (2)  classified  as  defective  is  (in  an  obvious  notation) 


j*jBinominl  0\,1\)  Y  multivariate  llypergeometri c  fn,N,N) 


where  Y  denotes  summation  subject  to  r,+  ...  +  r.  =  r  and  r 

L r  1  k  lrrr2 . 


is  the  multinomial  coefficient  r ! /  II  r.!  . 

In  particular,  with  p  =  N  [  N.p.  (again,  denoting  the  probability  that 

j=]  3  3 

an  individual  chosen  at  random  is  classified  as  "defective" 


h(Z)  =  up 

var(Z)=  np(l-p)  -  K^p.-p)2 

(Note  that  var(Z)  is  usually  less  than  the  "binomial"  value  np(l-p).) 


(0.1) 

(6.2) 


2.  GRADING 


Formulae  (4)-(6)  can  he  regarded  as  generalizations  of  (!)-(.>).  Further 
useful  generalization  is  obtained  by  considering  judgement  not  to  be  restricted 
to  "defective"  or  "non-defective"  but  to  include  assignment  to  one  of  several 
categories  -  as  would  be  the  case,  if  product  were  being  graded  in  terms  of 
quality  and/or  size  with  regard  to  marketing.  (One  aspect  of  this  situation 
is  the  multivariate  topic  of  disimminent  analysis,  but  here  we  are  concerned 
with  the  consequences  rather  than  the  methods  of  assignment.) 

Ke  will  analyze  a  situation  in  which  the  aim  of  judgement  is  assign 
an  individual  to  one  of  s  classes  C  ,  ....  C^.  We  denote  the  probnbilith 
that  an  individual,  who  really  belongs  to  C.  will  be  assigned  to  C.  by 
P.j ,  with,  of  course  ^  ^  1 V  ^  =  1.  Still  further  generalization  is  possible  by 

introducing  stratification  within  each  class.  This  leads  to  straightforward, 
but  notationally  complicated  elaboration,  and  will  not  be  pursued  here. 

Let  Yj  denote  t lie  number  of  individuals  belonging  to  Cj ,  in  a  random 

sample  of  size  n  for  a  lot  of  size  N  containing  N.  individuals  in  ( i  =  1 . 

pi  s 

s;  /^_jN^=N);  and  let  Z.  denote  the  number,  among  these  Y.,  assigned  to 


V 
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(i=l,  ....  s).  The  Yhas  a  multivariate  hypergeometric  distribution  with 
parameters  (n;N;N),  so  that 

P(Y)  =  /  n 


3=1 


In.' 

3 

Y. 

3 

/(^  (  I  Y.-n) 

j  =  l  J 


(7) 


Also,  given  Y,  Z.  =  (Z. . ,  ....  Z  .)  has  a  multinomial  distribution  with  parameters 

ij  sj 

(Y.,P.)  where  P.  =  (P..,  ....  P  .),  and 
J  ~J  ~J  ij  sjJ> 


s  Z.  . 


P  ( Z .  i  Y )  =  Y.!  T?  {P.  J/Z..n  (V  Z..*Y.) 

V  1  v  1  1  11  v  ^  -.1  T  ^ 


J  j=1  U  ij 


i=i 


i.l  3 


(8) 


and  Z^,  ....  Z^_  are  mutually  independent.  It  follows  that  the  joint  distribu¬ 
tion  of  all  the  Z’s,  Z  =  (Z^ ,  ....  !  )  is 


fy)  -IS  J  1  )  S  Z  .  • 

P(Z)  =  P(zjYjPfY)  =!‘\  n  \  N.1-1  n  (P  ’3/z..« 


'  j=i 


1  i=i  1J  JJ 


T9J 


where  Z.  .  =  J  Z.  .  .  (.Note  that  i  lie  V.,  are  ••  /  <  iw  'net!  by  tlu  Z.  .  ,  .1 
1J  1.1  ,i's  '  i.i ’s 

Pormally  wc  can  write  the  dist ribut ion  of  Z  as 

S  A 

Multinomial  (Y..P(|Y  Multivariate  llypergeometric (n ,N,N)  (10) 


and  it  might  be  called  a  "multivariate  hyporgeometric-multinomial  distribution. 
The  joint  factorial  moments  of  Z  can  he  obtained  from 


f 

s 

(  r .  . ) 
s  n 

r  ~ 

s 

s  (rii) 

■ 

s  (r..)s  r  ' 

E 

11 

n  z. . 

=  K 

V. 

n 

n  z. .  J  | v 

=  i 

e  n  y.  13  n  r.1- 

j-i  1J  , 

Y 

U=1 

i*i  J 

•  i  j  .  ,  n 
3=1  i=l  '  J 

fr..)  s  /  'r. .)s  r.  \ 

=  T  n  {n.  j  n  r.'-'j 
N(r--)  j=1  <  i=1  n  ' 


(in 


s  s 


where  r. .  =  J  r..;  r..  =  J  J  r.  . 

n  ;  l*  L-  *  ■» 


L=1  J 


i=l  j=l 


In  particular,  with  i  i  ijj  ^  j' 


E(Z. .)  =  nN  N.P. . 
13  3  13 


E(Z.  .Z.  ,  .)  =  n(2:)N?2)!>.  .P.  ,  ./Nf2^ 

13  i  '3  3  3  3  1  ’3 

E(Z. . Z.  . . )=  n(2)N.N. ,P. .P. , . ,/N(2) 

13  1  3  331313 


(12.1) 


(12.2) 


02.3) 


whence 


cov  (Z.^Z.,.)  =  -nN.{(N-n)N.+N(n-l)}  P.  P.  ,  .N'2  (N-l ) _1  (13.1) 


cov  (Z..,Z.,.,)  =  -nN.N. , (N-n)P. .  P.  . ,N”2(N-1)_1 
1 3  i  '  j  '  3  3  1313 


(13.2) 


Using  these  results  we  can  find  the  covariance  between  Z..  =  £  Z. . 

s  1  j=l  ” 

and  Z...  =  T  Z. , .  -  that  is  the  total  number  of  individuals  assigned  to 
1  i=l  1  J 

C.  ,  C_.  ,  respectively, 

cov  (Z...Z.,.)  =  -n{(N-n)P.P. ,  +  (n-l) (P.P. ,) }  (U 


where 


—  -1  r  -  -1  r 

P.  =  N  y  N.P. . ;  (P.P. ,)  =  N  y  N.P. .p. , . 

1  •  j  3  13  l  1  .  j  1  11  1  ’  1 


(i\  is  the  probability  that  an  individual  chosen  at  random  is  assigned  to  C.  , 

( P i P i , )  is  the  probability  that  an  individual  chosen  at  random  would  be  assigned 

to  Cj  on  one  judgement  and  C. ,  on  another.).  The  marginal  distribution  of 

2..  is  like  (4)  with  k  replaced  by  s  nml  u.  in-  I'...,  so  we  can  use  the 
1  v  ‘  >  •  1 .1 

corresponding  formulae  for  moments. 

3.  CROUP  SCREEN  1 NC 


Further  interesting  distributions  arise  in  connection  with  "group 
screening"  (Dorfman  (1943)),  in  which  groups  of  units  can  be  tested  for  the 
existence  of  one  or  more  defective  units  among  them.  This  can  be  practicable, 
for  example,  when  testing  liquids  for  presence  of  contaminants,  and  is  then 
suggested  as  a  possible  way  of  reducing  the  average  total  amount  of  testing. 
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Suppose  that  material  l'rom  units  is  mixed  and  tested  lor  presence  of 
"defective"  material.  If  a  negative  result  ("no  defectives")  is  obtained, 
no  further  action  is  taken,  but  if  there  is  a  positive  result,  each  unit 
is  tested  separately. 

Let  Pj ,  pj  denote  the  probabilities  of  obtaining  correct  or  incorrect 
positive  results,  respectively,  at  the  first  test.  As  before,  p,  p’  denote 
the  probabilities  of  correct  or  incorrect  positive  results,  respectively, 
when  units  are  tested  individually;  0,  N  denote  the  number  of  defective  units 
and  the  total  number  of  units  in  t he  population  respectively,  and  Y  denotes 
the  actual  number  of  defective  units  among  the  n^  tested. 

The  overall  probability  of  obtaining  a  positive  result  on  the  first 
test  is 

{l-P.in^Pj  ♦  I’Jiypj  -  p1  -  (Pj-pppJnj)  (l.V) 

where  P^n^)  =  (N-IV)  ^nl  is  the  probability  that  the  sample  contains 

no  defective  items. 

The  distribution  of  Z  can  be  represented  as 

A 

(Binomial  (IVY,p)  *B  i  nomi  a  1  (h'(n  -Y)  ,p  * ) )  V  Hyper  geometric  (iij.D.N)  (16) 

where  HI  is  an  indicator  variable,  defined  by 

I )  if  the  first  test  gives  a  positive  result 
W  "  { 

h)  otherwise 


Denote  that 


P(K=1)|Y>0)  =  P]  ;  P(W=1) | Y=0)  =  Pj 

An  explicit  formula  for  P(Z=z)  is 


07) 


r 


P(Z-z)  =  t, 


W,  Y 


'  ,lvY.  x,,  WY-x 


r  ,  w  \ .  x  , , 

L  (  x )P  (1-P) 


Av f rr j  -Y) 


z-x 


-  v  W(n.-Y)-z+x 

p'  d-p') 


where  ^  denotes  summation  over  (0,z-W(n^ -Y))  s  x  .  min  (IVY ,  -) . 


We  have 


E(Z(T)  |W,Y)  =  l  KjCKYl^^fKCnj-V)  ifr"l)p1p'r'1 

i=0  l1' 


Noting  (17)  and 


l)(n1-Y)(r"l)  |Y>0)l’(V  0)  =  ] 


’  (N-D) (r'1  Vn^1-1  for  i>0 


for  i  =  0 


(18.1) 


lnJr)  |  (K-D)  {N(r)  >_1  -lynj) ) 


0  for  i  -»  0 


E(Y^  (rij - Yj  (r'1^  |Y=0)P(Y=0) 


n j  r  ^  rQ  («i )  for  i  0 


(18.2) 


we  find 


w(t'<  -  "!r)  (pf^r1  v 


In  particular 


(19.1) 


E(Z)  =  n^PjP  -  (|»1-pj>i’,l’0(n] 


(19.2) 


=  nj  (I’jP  - 


wlu-rc 

and 


'(1)  ■  <>vr;"’  l"j! 


?  _ 


var (2)  = 


n."  /  »(p“-  p  “  )+(.v-n)i£,_-  p  “ )  i, 

-  M  .  . .  -  ’  “  r 


I  -  M  <— 

n  i  Lpiu’ 


N  (  N  -  1  ) 


Wi 


+  nl  fl,lP  '  '’(l)'”’  '  Vm1'  '  W' 


1 IP..V' 


If  there  is  no  possibility  of  a  "false  po  iti\e"  (so  that  |> •  p •  =  0)  wc-  find 


H  f  72)  =  lijPjpD/N 

n^pL'  P,pP 

w;>  -TT-d-V) 


n  Cn  - 1 ) pl> 

- {N’CD- 1 )  -  (N- 1  )DP  ) 

N-(N-J)  J 


I 


->)- 

In  general,  it  is  to  be  expected  (and  hoped)  that  p^  >  pj  just  as 
p  >  p',  since  we  would  expect  (hope)  that  the  probability  of  correct  decision 
would  exceed  that  of  incorrect  decision.  It  may  well  happen  that  p^  <  p 
since  detection  of  a  defective  may  be  more  difficult  with  the  mixture  of 
material  from  separate  units.  More  complicated  distributions  will  be  obtained 
if  it  is  supposed  that  p^  depends  on  the  value  of  Y  (  the  number  of  defective 

units).  It  does  not  seem  unreasonable  to  suppose  that  might  increase  with  Y. 

The  effectiveness  of  the  screening  procedure  is  measured  by  the  three 
quantities 

(i)  Probability  of  correct  classification  for  defectives  =  p^p  (20.1) 

(ii)  Probability  of  correct  classification  for  nondefective  = 

pJfalHi-p'p’)  ♦  (l-Po(n1))(l-P1P')  =  MPi-Pfi))?’ 

where  *  *  *  (n.-l)  (n.-l) 

P(D  =  (prpi)po(ni);  and  W  =  /cn-d  (20.2) 

is  the  probability  that  the  sample  contains  no  defectives,  given  that  one 
member  of  the  sample  is  nondefective. 

(iii)  Expected  number  of  tests  =  1  +  n ^  f Pq ^ ) P{  +  ^ (n] ) ]  (20.3 

■  1  *  "i (pi  -  pcn1 

of  course,  the  larger  the  values  of  (i)  and  (ii),  and  the  smaller  the  value 
of  (iii),  the  better. 

From  (20.1)  it  is  clear  that  the  probability  of  correct  classification  of 
a  defective  is  decreased  by  the  screening  process,  (since  p^p  <  p) .  The  value 
of  screening  must  therefore  come  from  increased  correct  classification  of 
nondefectives  or  reduction  in  the  expected  number  of  tests.  Table  1  contains 
some  relevant  numerical  information.  In  the  absence  of  screening,  Hj  tests 
would  be  necessary,  so 


L-U+n^Pj-P-j.  }/n,  =  1-n j1  -p^ 


indicates  the  proportionate  saving  from  screening,  and  this  is  given  (as  a 

percentage)  in  the  last  column  of  Table  1.  We  have  taken  N  =  °°,  so  that 

to  =  D/N  is  to  be  interpreted  as  proportion  defective,  because  the  values  do 

* 

not  depend  greatly  on  N.  As  N  decreases,  P^  and  P  decrease  so  that 
both  the  proportionate  saving  in  number  of  tests  from  screening  and  the 
probability  of  correct  identification  of  nondefectives  decrease,  (though 
not  substantially,  unless  N  is  quite  small).  It  is  to  be  noted  that  screening 
increases  the  probability  of  correct  identification  of  nondefectives. 


4.  HIERARCHAL  GROUP  SCREENING 


Sometimes  additional  saving  in  the  expected  number  of  tests,  and  improved 
accuracy  in  classification,  can  be  attained  by  using  two-or  more-stage 
screening  -  that  is,  hicrarchal  screening.  For  simplicity,  we  will  consider 
two-stage  procedures,  with  a  first-stage  sample  of  size  n^  =  hn^.  (Generaliza¬ 
tion  to  more  than  two  stages  follows  similar  lines).  If  a  positive  results  is 
obtained  for  the  combined  sample,  it  is  split  into  h  subsamples,  each  of  size 
n.  and  each  is  then  treated  as  in  Section  3.  Letting  p_,,  pi,  denote  the 
respective  probabilities  of  correct  and  incorrect  positive  results  when  testing 
each  of  the  second  stage  (size  subsamples,  the  three  quantities  measuring 
effectiveness  are: 

(i)  Probability  of  correct  classification  for  defectives  =  p^p^p  (.71.1) 

(ii)  Probability  of  correct  classification  for  nondcfectivcs 

=  P(*(nj)  (l-Pjp2P)  +  0-I,*(n1))(l-p1)  +  (P*(n2)-P*(n1))p1(l-p2p') 


+  (l-P0(n2))pJ(l-pJp’) 


I  *  f  1  »  *  ft 

=  1-P2P2P  +  lW(Pi-P1)P2P  +  P0<'n2)Pl  fI12_P2)P 
=  l‘P1P2P  +  C>’JnP2  +  P(2)P1)P  (21.2) 

★  *  t 

where  P  =  r0r J»2^  (P2"P2) 

*  * 

(Note  that  P^fn^  -  T^Cn^)  is  the  probability*  that  a  random  sample  of  n^ , 
known  to  contain  at  least  one  nondefective,  also  contains  at  least  one 
defective,  but  a  randomly  chosen  subsnmple  of  size  n,  containing  at  least  one 
nondefective,  in  fact  contains  no  defectives. j 
(iii)  Pxpected  number  of  tests 

=  1  +  h(Pj-P(n)  +  n1[P0(nJ)p]p2  +  (J-l’0(n2))Plp2  +  fp0rn2)-p0(ni))plp21 
=  1  4  h(prP(1))  4  n](P]p2  -  P(n?2  -  P(2)Pl)  (21.3) 

f 

where  P  ^  =  P0(nnl (p,-p,l . 

Th e  proportional  reduction  in  expected  number  of  tests  is 

1  -  Pjr_,  -  "21  lpJ  '‘\l  )  1  +  P(1)P2  +  p(2)pl  C;i  .4  1 

Values  given  by  (21.21  and  (21.1)  are  shown  in  Table  2. 

From  (21.1)  we  see  that  the  probability  of  correct  classification  of  ;. 
defective  is  decreased  more  than  for  simple  screening  (cf.  (20.1)).  As  some 
compensation,  the  probability  of  correct  classification  of  nondefectives  is 
higher. 

As  in  the  case  of  simple  screening  the  advantages  of  a  screening  procedure 

are  greater  when  the  proportion  fu)  of  defectives  is  smaller.  The  effect  of 

*  * 

finite  lot  size  (N)  is  to  decrease  1’.^,  P^,  P(]s,  P^ .  From  equations  (21) 
it  can  be  seen  that  this  will 

(il  not  affect  the  probability  of  correct  classification  of  defectives 
(ii)  decrease  the  probability  of  correct  classification  of  nondefectives 
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(iii)  decrease  the  expected  number  of  tests. 

Some  analysis  of  the  distribution  of  the  total  number  (Z)  of  items 
classified  as  "defectives"  by  this  hierarchal  screening  procedure  is  given 
in  the  Appendix.  Here  we  just  give  the  expected  value 


,-l 


E(Z)  =  n1[p1p2(p-p  )DN~  +  (P]p,  -  p(2)1'i'l>(]  jl12')p  1 


C2j 


AC KNOW I  i.DriIMF.XT 

Samuel  Kotz's  work  was  supported  by  the  II. s.  Office  of  Naval  Research 
under  Contract  N0O014-81 -K-0S01 .  Norman  I.  Johnson's  work  was  supported  by 
the  National  Science  Foundation  under  Cram  MCS-8021704. 


PROBABILITY  OF  CORR!  CT  1YPLCTED  PERCENT 

CLASSIFICATION  OF  NONDCFF.CTIVE  REDUCTION  IN  TESTS 


PROBABILITY  OF  CORRI.CT  FXPFCTliD  PtikCLXT 


H 

a 

K- 

a* 

o 

*— < 

•—  — 

m 

K-* 

3 

►— 

Table  2:  Two-Stage  lliorarchal  Scrcenin 


PROBABILITY  OF  CORRECT  IXPLCILU  PLRCLXT 
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ln  order  to  study  the  distribution  of  the  total  number  (Z)  of  items 
declared  "defective"  we  introduce  the  auxiliary  indicator  variables: 


w  _  J  1  if  items  in  the  j-rh  subsample  are  tested  individually 
j  \  0  otherwise 

(Al) 


for  j  =  1,  h.  Conditional  on  W  ,  . . . ,  and  the  actual  numbers 

Yj,  ....  of  defectives  in  the  corresponding  subsamples,  Z  is  distributed  as 

h  h  , 

Binomial  (  \  W.Y.,p)  *  Binomial  N. (n^-Y . ) ,p  )  (A2) 

j=l  3  3  j«l  3  ~  3 

The  conditional  r-th  factorial  moment  of  Z  is 


C(Z(r)|W,Y)  = 
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We  have 


PCY  =  0)  =  P0(n2);  P(Y.>0)  =  1-P0(nj 


P(Yj  =  0,  Y  •  =  0)  =  PQ(2n2) 


P(Y.  =  0,  Y  •  >  0)  =  P(Y  ■  >  0,  Y.  =  0)  =  p0(n2)-p0(2n2) 


PCYj  >  0,  Y  •  >  0)  =  l-2P0(n2)  +  p0(2n2) 

P(Yj  Y  *  °*  jji  >  «  "  V2^  '  W 

h  h 

P(Y  =  0  =  l  Y  )  =  P(Y  =  Y  •  =  0  =  l  Y  )  =  Pn(n  j 
3  i=l  J  3  i=l 


whence 


E0v'“)  =  PlP2  -  p(2)Pl  -  P(1)P2 

E(wJ  w“,)  =  Plp2  -  2P]p2P(2)  ♦  P](p2-p2)2P0(2n2) 


O) 


(A6.1) 

(A6.2) 

(A6.3) 

(A6.4) 

(A6.5) 

(A6 . 6) 

(A7.1) 

CA7.2) 


Also,  from  (A4)  and  (AS),  with  6,  fi '  >  0 


e(wJy?)  =  Pip2ecyJ) 

E(WaWa,,YeY(5')  =  p  p2E(YeYe') 
J  3  J  J  1  l1  2  J  J 


(A8.1) 

(A8.2) 


ECvAi^yJ)  =  plP^{l-P0(n2)|-  E CY J | Y .  ,  >  0)  ♦  lYj  »“0)  (A8.3) 


v8 


B.v 
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The  joint  distribution  function  of  V  is 
PCV)  = 


CA9.) 


D 

N-D  | 

1  /  N  1 

h 

h 

/ 1 

vi . v D  -  iYi 

^  ) 

VY1»  •••>  n2~Yh ’  N-D-nl+ 

1  J 

n 

n2  ’  •  •  •  *  n2 >  N-n j  j 

and  the  r-th  joint  factorial  moment  of  V  is 


n  y. frj1'! 

U-i  -  J 


(  h  , .  , 

TJ  fl  ‘  *  i  / 


j  =  l 


11  J 


I)(r)/Nfrl 


where  r  =  L  r..  In  particular  for  all  j  ^  i 
i  =  l  1 


(A10) 


F  (Y  )  =  nJ'/.Y 

t'Al  1.1} 

H(Yj‘)=  n^P  in  1  ,  ■  / \ *  "  1 

(AH  .2) 

1  !  M  ("M 
' . Y .  * )  =  n ~ U  “  ,/NUJ 

J  J  - 

(All. 3) 

=  0)  =  n,D/(N-n,l  ;  I  (V.  X  • 

°>  =  n2D{ 3. -rt2N~ J  _p(i  fn^) .  1  -P( 

,  (n-,)  ;_1 

*  (Nnj'1 

(All .4) 

J  J 


Applying  formulae  (A3)  -  (Ml'  wo  get  alter  some  algebraic  rearrangement  the 
formula  (22)  for  F.(I).  In  order  to  calculate  the  variance  of  Z,  we  have  to 

calculate 


E(Z^2))  =  H[p‘(  J  W.l.)(  J  bjlj-ll  +  2pp  (  7  IV  -  Y .  )  (  l  w.(n,-Y. 

j  =  ]  -  -  j  =  ]  -  i = 1  i  =  ]  J  “  1 


,  h 


-VJ) 


•  f  f  h 

M>  t  7]1'  j  ,HJ  '  i !  1 1  _L  IV.  (ii,-Y.)-  1)J 


Cp-p  )  E[C  l  IV.  V  )  ]  -  (p  - 
3-1  3  3 


P  U:<  L  h;Y.) 
i  =1  J  J 


,  •  h  h  t2 

+  2n.p  CP-P  )E[(  I  W.V .)(  I  W.J|  -  n,p  E(Vtt.) 

j=l  J  J  j=l  J  J 

2 

♦  n^p'  EfC^V.r] 


Cp-p,)2[hF(lV.Y.2)  +  h(h-l)E(W  .l\\  'Y.Yj  ')  1  -  p2-p,2)hE(W.Y  ) 

it  ,2 

+  2n.,p  Cp-p  )  £hE (W .  Y  j  )  +  h fh- 1 lEOV .hZ  'Y^ ) ]  -  n _,p  hE(lV^) 


+  n;p  [hF.CK.1  +  h(h-l)i:(W  .W  ' )} 


where  E(l\Z)  and  HflV.lV.')  are  given  by  (A-  i,  a; 


E°'jYj)  =  njPlp2nN" 

E(K.Y.2)  =  n 7p^p^D{ (ji.-l  1  D+N-n,  'N* 1  (N  - 1  !  "a 
H(W.IV.»Y.)  =  n2r1p2|l'N"lp.’-(’N‘M2l"1|’(:i  J 


The  resulting  formulas  for  EtZ  ”  j  and 


var(Z)  =  E(I  ' ' )  +  E(Z)  -  1 1 : ( Z 1 


are  complicated,  but  numerical  calculation  is  straightforward.  If 

1 

p;  =  0  (so  that  there  ire  no  false  posit i\cs),  then 


E(Zl“J)  =  p“[hl:(lV  Y.'t  +  ii  ( h  1  )  i:  ( l\  l\  .  M  i  .  '  >  -  hlOV.Y.  j| 
1  1  I  J  i  j  .1  j  j  i  1 


njPjP2p-b  11(1'  1  ) 


(Ii- 11  ip  1  1  -n,.\'  )P 


and  also 


*j<m  & 


.'V.*  ■ 


NOTATION  SUMARY 


Single  Sampling 

p  Probability  that  a  defective  item  is  classified  as  defective 

p*  Probability  that  a  nondefective  item  is  classified  as  defective 

p  =  N  ^Dp  +  N-D)p'}  =  Probability  that  an  individual  chosen  at 

random  is  classified  as  defective. 


Hierarchal 

Sampling 

n . 

J 

Pi 

Number  of  items  tested  en  bloc  at  j-th 

stage 

Probability  that  a  group  containing  at 
classified  as  defective,  at  j-th  stage, 

least  one  defective  is 

V 

Probability  that  a  group  containing  no 
as  defective,  at  j-th  stage 

defectives  is  classified 

vv = 

Probability  that  a  random  sample  of  n. 
items  (i  =  1,  2).  •' 

items  contains  no  defective 

Vnj>  ' 

Probability  that  a  random  sample  of  n. 
items,  given  that  it  contains  at  leas^ 

items  contains  no  defective 
one  nondefective. 

P(j)  = 

ciypjJ'yv 

* 

p(j)  ’ 

cpj-p^WV 
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